Kondo effect in a double quantum-dot molecule under the effect of an electric and 

magnetic field 
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Electron tunneling through a double quantum dot molecule, in the Kondo regime, under the effect 
of a magnetic field and an applied voltage, is studied. This system possesses a complex response 
to the applied fields characterized by a tristable solution for the conductance. The different nature 
of the solutions are studied in and out thermodynamical equilibrium, ft is shown that the interdot 
coupling and the fields can be used to control the region of multistability. The mean-field slave-boson 
formalism is used to obtain the solution of the problem. 
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The electron tunneling through a double-quantum-dot 
(DQD) molecule in the Kondo regime has received much 
attention in the last yearsiSi 3 ^^^*^!!^ Recently Jeong et 
alii reported the observation of a coherent Kondo effect 
in a double quantum dot, which is a direct experimental 
realization of a two-impurity Anderson model^^*^*^ 
The strength of the coupling between the QD's deter- 
mines the character of the molecular electronic states 
and the transport properties of the artificial molecule. 
In the tunneling regime, the states are extended across 
the entire system and can be constructed from a coher- 
ent state based on the bonding or anti-bonding levels of 
the quantum dots. These many-body molecular states 
are observed in a transport experiment as a differential 
conductance of two peaks, a consequence of the splitting 
of the Kondo resonance. 

Recently Aono and Eto studied theoretically the prob- 
lem of a DQD in thermodynamical equilibrium, under 
the effect of a magnetic field oriented parallel to the cur- 
rent directioni&They calculated the magnetoconductance 
and found a suddenly drop of the conductance as a func- 
tion of the field, when t c > T, where t c is the interdot 
interaction and T the resonant width of the dot levels due 
to coupling among them and the contacts. Some aspects 
of the non-equilibrium transport properties of a DQD in 
the Kondo regime, under the effect of an applied voltage, 
has recently been studied^. When the inter-dot coupling 
is greater than the level broadening and there is an exter- 
nal voltage applied to the system, the conductance shows 
a tristable behavior. 

In this work we analyze the physics derived from the 
complex behavior that characterizes the electron tunnel- 
ing through a DQD molecule in the Kondo regime under 
the simultaneous effect of a magnetic and an electric field. 
We show that the discontinuity of the conductance found 



by Aono and Etoj£ increasing the magnetic field when the 
system is in equilibrium, is a result of a tristable behav- 
ior of the magnetoconductance. This results is similar to 
the behavior the system has when it is out of equilibrium, 
due to the application of an external voltage^ 

We consider a DQD connected in series attached to two 
leads. The physics associated to this artificial molecule 
electrically connected can be readily understood in terms 
of two impurities Anderson Hamiltonian where the im- 
purities located in the site and I are the quantum dots. 
The leads are represented by a one dimensional tight- 
binding Hamiltonian connecting the dots to two particle 
reservoirs characterized by Fermi levels fir, and fj,R re- 
spectively. The left (right) dot is connected to the left 
(right) lead by the hopping Vl(r) and the dots are con- 
nected among themselves by the interdot tunneling cou- 
pling t c . Each dot a has a single energy s a (a = 0,1). 
Under an applied magnetic field the energy levels of the 
dots are split by the Zeeman effect into e a ± agfisB. We 
adopt the two-fold degenerate Anderson Hamiltonian in 
the limit U — > oo, that is diagonalized using the mean- 
field slave-boson formalism>i2iii*i2i The double occupancy 
in each dot is forbidden and the inter-dot Coulomb in- 
teraction neglected. We introduce the slave-boson oper- 
ator &t that creates an empty state and a fermion op- 
erator f aa that annihilates a single occupied state with 
spin a. To eliminate the possibility of double occupancy 
we impose the constraint Q a = Y^ a f\ a faa + b^ a b a = 1. 
The annihilation operator of an electron in the dot a 
is c Q(T = w a f aa . In the mean field approximation the 
bosonic operators b\ and b a are replaced by their ex- 
pectation values, (b a ) = 6 ce v / 2 = = b\^/2. Hence 
the Hamiltonian of the DQD molecule connected to leads 
plus constraints is written as, 
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a— 0,1, cr (7 a 

+ to Y^flha + flM + E A « ffifia - 1)- (1) 



a=0,l 



where e a = e Q + X a , V L = V L b , V R = V R bi, t c = t c b b ll 
the X a are Lagrangian multipliers, which guarantee the 
constraint conditions on the Q a . 

The first term on the right-hand side of Eq. 1 represents 
the electrons in the left and right leads; 

Hl ead ^ £ i n iT+ t E 4a C 3<y ( 2 ) 

cr,i^0,l <ij^0,l><T 

where the operator cj CT creates an electron in the site i 
with spin a, Si is the site energy and t is first-neighbor 
hopping in the leads. Here Hi ea d is a free-particle Hamil- 
tonian with eigenfunctions of the Bloch type, 



The stationary state of the complete Hamiltonian H 
with energy e k can be written as 

IVW = ^a k a \j,a), (5) 
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where the coefficients a k a are given by, 



a% =< j,o\i>k„ > ■ (6) 

We obtain the following eigenvalue equations for the 
Wannier amplitudes a k a 



\k,a)=Y,e ikja \j,a), (3) 
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where \k,a) is the momentum eigenstate with spin a and 
\j) is a Wannier state localized at site j of spin a. The dis- 
persion relation associated with these Bloch states reads 

e = 2tcos(kd). (4) 

The model has an energy band, extending from — 2t to 
+2t, with the first Brillouin zone expanding the k interval 
[— n/d, Tr/d]. 



™ 3 ,o = (j,v\H\ip k a) 

ea\ a = e ja % a + t(a]_ l cy + a k +la ) (j £ -1, 0, 1, 2), 

£CL -l(2),a = £ -l(2) a -2(2),<r + V L(R)%(1) ><7 + to -2(3),<r' 

£a a(i),a = (£o(i)+vgVBB)a* {1)a + V L{R) a k _ 1{2)a 

+ * c <o),<r> (7) 
where a k a is the probability amplitude to find the elec- 
tron in site j and state k with spin a. 

The four parameters (bo, b\, Ao, Ai) are determined 
minimizing the expectation value of the Hamiltonian (1). 
They satisfy the set of equations: 
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(8) 



k.a 



The sum over the wave vector k cover all the occupied 
states. The resulting equations are nonlinear due to the 
renormalization of the localized levels in the dots, the 
interdot coupling tunneling and the coupling tunneling 
between the quantum dots and the leads. 

In order to study the solutions of equations (7) and (8) 
we assume a plane wave incident from one lead with an 



amplitude / and with a partial reflection amplitude R. 
The waveform at the other lead is a simple plane wave 
with intensity given by the transmission amplitude T. 
Taking this to be the solution of the system at infinite, 
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for positive fc, we can write 

J _ T+Jkdj +R + e -ikdj ^ ■ <Q 



(9) 



and those with negative fc 



= I-e 



-ikdj 



R- 



Jkdj 



"jcr o J 



J>1, 

i < o. 



(10) 



These functions describe free electrons approaching the 
scattering potential from j = — oo and j — +oo, respec- 
tively, with transmission and reflections amplitudes T+ 
and .R+ for k > and T" and R~ for fc < 0. 

The solution of equations (7) can be obtained through 
an adequate iteration from right to left for k > 0, and 
from left to right for k < 0. For a given transmitted 
amplitude, the associated reflected and incident ampli- 
tudes may be determined by matching the iterated func- 
tion to the proper plane wave at the left lead for fc > 
and at the right lead for fc < 0. The transmission am- 
plitude t a — T+/I+ and t a = T-jl^ are obtained using 
an iterative procedure. The total transmission proba- 
bility with a spin a and energy e — 2t cos(fcd) given by 
T(e) = ^(e)! 2 . Then the linear conductance at 

zero temperature is written as, 



2<= 2 
G=^T. 
h 



(11) 



The equations (7) and (8) are nonlinear and require a 
self-consistent solution that is obtained using a conjugate 
gradient algorithm. In this work, we are interested in an- 
alyzing all the stationary solutions for the conductance. 
As G(£?)and B(G) are multivalued functions, the multi- 
ple solutions for the quantities 6o, &i, Ao, Ai are obtained 
by fixing the magnetic field B and obtaining the value 
of the conductance G that corresponds to the solution 
of the equations of motion of the system. This method 
permits us to obtain the complete G — B curve. 

We study first a model which consists of two leads 
equally connected (Vl = Vr = Vq) to both sides of a 
DQD dots in equilibrium with chemical potential /i^ = 
fi R = 0, t - 30 T , V = 5.48 T , e = e 1 = -3.5 T 
(Kondo regime with T K » 10~ 3 T with T = 2 7 rV 2 p(0)). 

Fig 1 shows the linear conductance versus the applied 
magnetic field for different values of t c with t c < To. For 
small values of t c , t c <C Tq our results agrees with those 
obtained by Aono and Eto. The magnetoconductance 
is negative in this region. However for greater values of 
t c , t c ~ Tq another behavior of the magnetoconductance 
starts. Within a range of values of the magnetic field 
there is a tristability behavior with three stationary so- 
lutions of the conductance for each value of the magnetic 
field. Fig 2 shows the results for t c > Tq. As in the re- 
sults found by Aono and Eto, at a critical value B c ^ there 
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FIG. 1: G-B curve for a)t c = 0.2r , b)t c = 0.6r , c)t c = 
0.8r and d) tc = T 



would be a discontinuous reduction of the conductance 
when B is increased and other B c i(B c ^ > B c |) where 
the conductance would rise discontinuously, when B is 
reduced. Using this procedure, the interior of the trista- 
bility remains inaccessible to the experiment. However, 
the internal shape could eventually be reached, as it is 
for the case of double barrier heterostructrures, if it were 
possible to employ a load line in the experimental mea- 
surement to go along all the points of the G — B curve. 
From this point of view the electrons goes along a non- 
linear potential, that is magnetic field dependent. In this 
case the magnetic field plays the role that the external 
electric potential plays in the case of a standard double 
barrier heterostructure, where the current has as well a 
multistable behavior^ 
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FIG. 2: Linear conductance G versus the Zeeman splitting 
gii B B for a)t c = 1.0r , b)t c = 1-2T , c)t c = 1.6T and d) 

tc = i.8r 

Increasing the magnetic field the tristability probably 
arises from a many-body ground state that changes its 
character due to a energy crossing. It is a singlet of the 
Kondo type when there is no applied magnetic field, while 
for a critical field onwards the ground state results to be 
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the component of the triplet wave function, (two paral- 
lel spins one at each dot) with spin projection parallel 
to the magnetic field. In the parameter space this two 
ground states have their own basin of attraction consti- 
tuting metastable states when they loose the condition 
of been the ground state of the system. As soon as the 
metastability is lost increasing or decreasing the field, the 
system abruptly relaxes to the ground state producing a 
discontinuous behavior of the current^ 

It is very interesting to study these non-linear effect 
when the system is under the simultaneous action of an 
external magnetic and electric field. In this case the DQD 
is out of equilibrium as it operates within two chemical 
potentials \il = V/2 and (1r = —V/2 defined by the 
applied dc bias voltage V. Imposing a fixed value of the 
magnetic field, once the amplitudes a*- a are known, the 
current can be numerically obtained from, 
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FIG. 3: J-V characteristic curve for t c = To for various values 
of the Zeeman splitting g\i B B. &)gyL B B = 0.0 , b)g/j, B B = 
0.407V, c)gn B B = 0.467V and d) g/i B B = 0.477V 

We study first the problem when t c = T. For this value 
of t c , in the absence of an external magnetic field, the sys- 
tem is at the threshold of having a manifestation of the 
non-linearity produced by correlation. For a regime in 
which g^sB <C Tk, figure 3 shows that the magnetic 
field reduces the height of the peak and valley of the 
J — V characteristic curve. However, for larger magnetic 
fields, < Tk, a tristability behavior of the cur- 

rent is promoted. For greater values of t c (t c = 1.4r), 
where the system has already a highly non-linear behav- 
ior in the absence of magnetic field, the shape of the 
tristability is highly dependent on B. As shown in fig- 
ure 4, here as well, a remarkable increase of the peak 
to valley ratio is obtained in the regime < TV/. 

This ratio goes from about 20 at zero magnetic field to 
about 2000 at g^sB — 0.333Ta'- These phenomena are 
a result of the reduction and eventual elimination of the 



Kondo correlation when the magnetic field is of the order 
of the Kondo temperature. The Zeeman effect destroys 
the coherent behavior of the two QD's produced by the 
molecular Kondo coupling, the electrons localize at each 
dot occupying levels that are mismatched by the applied 
potential and consequently the current goes to zero. 
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FIG. 4: J-V characteristic curve for t c = 1.4To for various 
values of the magnetic field. a,)g/j, B B = , b)g/j, B B = 0.27V, 
c)g^ B B = O.STk and d) gf i B B = 0.337V 

In non-equilibrium the nonlinear behavior described 
above has similarities with the well-known bistabil- 
ity present in three dimensional resonant tunneling 
double barrier systems and multistabilities in doped 
superlatticesi2iiiiSiiSiiLi£ In three dimensional tunnel- 
ing double barrier systems, the Coulomb interaction pins 
the renormalized energy level at the well when the ap- 
plied voltage V is augmented maintaining the system 
at resonance. However, when V goes above a critical 
value, the abrupt leakage of the charge accumulated in 
the well takes the system out of resonance. For doped 
superlattices, at large carrier densities, the internal po- 
tential profile produced by the applied voltage breaks up 
into domains giving rise to a multistable behavior in the 
current ^i&ii I n both systems, the non- linear deforma- 
tion of the potential profile due to the electron-electron 
interaction is responsible for the phenomenon. Although 
in our case the Kondo effect is a strong many body effect 
of a different nature in comparison with the previous sys- 
tems, the multistability is as well a consequence of the 
non-linear potential profile deformations that the elec- 
trons suffer when the external potential is modified. 

In conclusion, we have investigated the magnetocon- 
ductance in a DQD molecule in the Kondo regime. We 
have shown that the magnetoconductance has a tristable 
behavior depending of the ratio of t c /T and the value 
of the magnetic field that controls the non-linearity of 
the system. This evidence was obtained using a mean- 
field approach that neglects fluctuations, appropriate for 
large spins as it is the first term of an l/N expansion, 
where N is the spin dimension. Some experiments have 
shownii that it is possible to design dots with a configu- 
ration that, due to the exchange interaction, the parallel 
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spin coupling following Hund's rule gives rise to a total 
dot spin S t > 1/2. For this cases our results are more 
reliable as fluctuations are less important. Although we 
believe that our results are robust against them, a study 
that goes beyond mean field is necessary to confirm our 



conclusions. 
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